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1 Introduction 



Recently, there is growing interest in difference versions of Painleve equations [T|]. A first ex- 
ample, the discrete Painleve I equation (dPI), which is a nonlinear non-autonomous ordinary 
difference equation of the form 

-±— = (3 (R n+1 + R n + i? n _0 + 7 , (1) 

tin 

appeared in the theory of matrix models for 2-dimensional quantum gravity, cf. There 
it occurs as lowest non-trivial term in the so-called 'string equation'. 

Other discrete Painleve equations were subsequently found. A discrete version of Painleve 
II was obtained in 0, starting from unitary rather than hermitean one- matrix models, and 
independently in [0 in connection with similarity reductions of integrable lattice equations. 
Next, using the new method of singularity confinement, ||, discrete analogues of Painleve 
III, IV and V, were found in 0. As the singularity confinement approach is based on what 
is still a conjecture (namely, the non-propagation of spontaneous singularities in integrable 
mappings), proofs of integrability of these equations had still to be provided. For dPI and 
dPII, the situation was clear from the beginning, because their very construction (either by 
the approach of orthogonal polynomials or by the similarity approach) led to the existence of 
isomonodromic deformation problems for these equations, cf. || respectively @. However, 
for the new discrete equations, dPIII-dPV, a priori no Lax pair was known. For dPIV and 



dPV the situation is still open, but in |10| a Lax pair for dPIII was presented (among new 
isomonodromy problems for dPI and dPII, and their variants). Interestingly enough, this 
Lax pair is a discrete isomonodromic deformation problem of a linear q- difference equation, 
rather than of a differential equation. It is to my knowledge the first time that g-difference 
equations arise in connection with Painleve equations. To study the asymptotics of such 
type of systems, it is necessary to go back to the historic works of the Birkhoff school on 
this subject, cf. e.g. [|11]]-|TJ|], that unfortunately seem to have fallen into oblivion since the 



1940's. Although it is relatively easy to find special solutions for the dPIII equation, cf. 
|T6|| , the analysis of the corresponding g-difference system seems quite complicated. As the 
theory of isomonodromic deformations for g-difference systems still needs to be developed, 
it would, be useful to have at ones disposal some slightly less complicated systems than the 
one for dPIII, in order to study the features of this kind of analysis. 

In this note, I would like to introduce for that purpose a new discrete Painleve type 
of system, namely a g-deformed version of the dPI equation ([]]). This is probably the 
simplest g-deformed transcendent, although it is naturally written as a third order difference 
equation rather than a second-order one like (|IJ). Nevertheless, the g-deformed discrete 
Painleve I equation, carries -like dPIII- an underlying monodromy problem of g-difference 
type. This system is a natural candidate to be investigated from an analytical point of view, 
namely by constructing asymptotic solutions and the connection coefficients between them 
and formulating an inverse (Riemann- Hilbert type) of problem along the lines laid out in 



T2|j . The isomonodromic deformation can then be used to 'solve' the corresponding Painleve 
transcendent. I will leave this program for a future publication, and in this note I will confine 
myself to establishing a connection between the g-deformed Painleve I and g-orthogonal 
polynomials, very much along the lines in which this connection was made in the one-matrix 
models. However, in order to do do this, we will notice that one needs to modify slightly 
the equation, namely by introducing one term higher in the string equation. Furthermore, 
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by making this connection we will be naturally led to some discrete-time analogue of the 
Volterra equation, the similarity reduction of which is the g-deformed Painleve I equation. 



2 g-Hermite Polynomials 



(/-Analysis, p~7|, 113] is very fashionable nowadays, mainly because of its use in quantum 



groups and the investigations on g-special functions in connection with their representation 



theory, cf. e.g. |0|-E2[]. In this note I would like to draw attention to a new class of g-special 
functions, namely g-difference deformed versions of Painleve transcendents, in particular a 
g-deformed version of the dPI equation (|I|) . It is well-known that the latter equation arises 
as the lowest non-trivial term in the string equation for one-matrix models via orthogonal 
polynomials. A special subcase, namely if a = (5 = in (0), gives rise to the Hermite 
polynomials. In seeking a g-deformed version it is natural, therefore, to start from the basic 
g-Hermite polynomials, defined e.g. in |17| . Let us first review these, mainly for the purpose 



of fixing the notations. 

First, let me recall the definition of the g-exponential functionF] 



00 r n / \ — l 



e x 

q — in, 



introducing also the g-numbers (n) q = ^— p • Sometimes it is useful to work with Andrew's 
notation, fT8|j , 



j=0 

in terms of which we have e.g. 



(ag n ; g) c 



e»-! = ((g-l^g)^ , (g<l). (3) 
Now we can introduce the basic g-Hermite polynomials by means of the generating function 



>3C 



-in 

S q (x,t) = efef. 2 = ^ T-^Hnix) , (4) 

Furthermore, introducing the g-differentiation and the g-dilatation operator 

gD x f{x) = — , q T x f{x) = f{qx) , 

(q — l)x 

and recalling that the g-exponential obeys the difference rules 

qD x &q , qPx ^q- 1 ^q -1 J 

one can derive for S q the following relations 

qD x S q = tSq , (5a) 
qD t Sq = (X+ (2) q at) q T~ X q T t Sq , (5b) 



1 In contrast with the original exponential function, the g-exponential function e x q has a finite radius of 
convergence R = |1 — ql^ 1 as \q\ < 1, and simple poles at x = (1 — q)~ 1 q~i for j = 0, 1, 2, . . . . 



3 



leading to the relations 

qD x H n = (n) ? -ff n _i 



H n+l =xq n q T- l H n + a(2) g (n) ? g n - 1 ^- 1 i? n _ 1 . 



(6a) 
(6b) 



The last equation can be converted into 

H n+ i = xH n + a(2) q -i(n) g q n H n _ 1 , (6c) 

whereas, from (|6a[) and (6b") one can also derive the second order g-difference equation 

q' n {n) q H n = x q D xq T~ l H n + a(2) q -i gD x P' 1 q D x H n , (6d) 

i.e. the g-analogue of the Hermite equation. Using the fact that qT~ x gD x = q q D xq T~ x , and 
that qD x = (2) q x q zD x 2, we can convert the last equation into 



—aq 2 x 2 

V 2 



F^tPsHn] + a(2) q -iq- n (n) q e-«fH n = 



(7) 



identifying a = (—a(q 1 + l) 2 ) 1 , which we take to be real positive. From (0) it is clear 
that the basic g-Hermite polynomials are orthogonal with respect to the Jackson integral 



f(x)d q x = (1 - q)c £ q n [f(q n c) + f(-q n c)] 

n=0 



1 - q 2 )a 



-1/2 



(0 < q < 1) 



2 

with weight function w g (x) = e~f% . Introducing the following representation for the q- 



Gamma function, [17, [L8|, 



r 9 (x) 



1/(1-9) 



i a q -id q t 



(g;g) a 
(?*;?)c 



:i-g) 



1-x 



in terms of which we can express the normalization constant, we have the orthogonality 
condition 

d a xeZ% H n (x)H m (x) = \ ^T q2 (-)^, m . (9) 



Finally, the dependence on the variable a (or, equivalently, on a) is significant. In fact, from 

t 2 



2D a S q 



(2W 95 



we derive 



q' 2 D a H n 



(n)q(w- l) g 
(2) 2 a 2 



n-2 • 



(10) 



3 A deformed Painleve I Equation 

Comparing the Lax pair for dPI, cf. e.g. 0, with the one for the g-Hermite polynomials, 
namely eqs. flSaf ) and (|6^), it is natural to pose the following g-linear system 

xP n {x) = P n +i(x) + i2 n P„_ 1 (a;) , (11a) 

qDxPnix) = A n P n „ X {x) + B n P n ^{ X ) , (lib) 
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the compatibility of which leads to the set relations 



A n+ i = qA n + 1 , 
B n+ i + R n A n _i = q (A n R n _i + B n ) , 
RnB n -i = qB n R n _3 . 



Eqs. (|12a|) and ( |12cj ) are readily solved as 

A n = (n), + ag n , S n = ^q~ n R n R n _ x R n _ 2 . 



(12a) 
(12b) 
(12c) 



(13) 



Inserting these expressions into ( |12b| ) we obtain the following third-order nonlinear non- 
autonomous ordinary difference equation 



f3q n (q l R n+ i - qRn-2) = q 



ji)q + aq r 

Rn 



(n-1) 



aq 



n-l 



Rn-l 



(14) 



which we will refer to as qPI. The connection is made clear by rewriting (113) in the form 



q 



ra-2 

7 + Pq~ 2 (Rn+l + Rn + Rn-l) + P(q' 2 ~ 1) £ Rj 



n 



+ aq n 



R„ 



(15) 



namely by performing one (formal) 'integration'. The infinite sum term in ( |15|) will clearly 
disappear in the limit q — > 1, in which case we immediately recover dPI , ([!]). 

In |23|] a g-deformed Painleve equation related to a special subcase of the discrete PII 
has been presented which shows some similarity to eq. (|14|). It is not clear, however, 
whether their equation can be reduced to flll]) by coalescence. Furthermore, the approach 
in p3| does not focuss on isomonodromic deformation problems, but rather starts from 



factorizing Schrodinger operators. What is interesting about (JT^), apart from being probably 
the simplest g-deformed transcendent that is now at our disposal, is that there is a connection 
with g-orthogonal polynomials, along the same lines as for the original dPI. However, in order 
to establish this connection one has to modify slightly eq. (|T4] ) as we shall see. 



4 (/-Orthogonal Polynomials 

We look for an explicit realization of the qPI equation in terms of g-orthogonal polynomials, 
in much the same way as the original discrete Painleve I equation (fl|) is solved by means 
of orthogonal polynomials. It is this solution that gives rise to a connection with discrete 
hermitian one-matrix models. 

For this purpose let us generalize the orthogonality conditions (Q) to 

dqXe^-2 e <j-4 Pn(x)Pm(x) = h n S nrn . (16) 

For the integration limits in the Jackson integral (11^) we take the smallest zero of the weight 
function w q (x) = e~_f e'-t , i-e. ±c, where c = min |((1 — g 2 )a) -1 ^ 2 , ((1 — g 4 )&) , 
(for fixed positive a, b, taking < q < 1) to ensure positivity of the measure. The polynomials 
are supposed to be normalized such that 

P n (x) = x n + ■ ■ ■ => <p x P n = (n) 9 P„_i + • • • . 
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Of course one may go further and introduce in ([16|) a weight factor depending on an arbitrary 
number of g-exponential factors, but I will refrain from doing so in the present note. 

From the orthogonality condition (fl6|) we can derive the isomonodromy problem ( |lla|) . 
with R n = h n /h n -i (as in the usual case). However, ( |llb| ) needs to be modified as follows 
from the following relation for the weight function 



<p x w q (x) 



{2) q ax + (4) q bx 3 + (q- \){2) q {A) q abx b w q {qx) 



which indicates that one has to push already to fifth order terms in order to take into account 
second and fourth order coupling constants a, b. This means that eq. ( |llb|) is not applicable 
here, and performing the calculation 



Jd q xw q (x) [( q D x P n )P m + P n ( q D x P m )] + (q-l)J d q xw q (x)x( q D x P n )( q D x P n 

= J d q xw q {qx){ q T x P n ){ q T x P m ) [{2) q ax + {A) q bx 3 + {q - l){2) q (A) q abx 5] 

we are led to a linear problem consisting of (|lla|) together with 

gD x P n (x) = A n P n ^(x) + B n P n ^(x) + C n P n _ 5 . 
Compatibility will now lead to 

A n+ i = qA n + 1 , 
B n+ i + R n A n _i = q (A n R n _i + B n ) , 
Cn+i + RnB n -i = q (B n R n s + C n ) , 
Rn,C n -i = qC n R n _5 . 

Eqs. (|19a|) and ( 19d ) again are readily solved as 



(17) 



A, 



(n) q + aq n , C n = 74 n R n R 



and from ( |19b| ) and (|19c| ) we get the system 
(n) q + aq n (n - l) q + aq n ~ x 



Q- 



Rn 



Rn-l 

B>„ — B. 



1 2 Rn+l — Rn- 



q 1 Rn+\B n+ i — qR n - 2 B 7 



n— 1 
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(19a) 
(19b) 
(19c) 
(19d) 

(20) 

(21a) 
(21b) 



in which B n = g~ n 7i?„_R n _ii? n _ 2 -Bn- Of course, one can derive from (|2l|) a closed sixth order 
difference equation by eliminating the B n . This is not so illuminating, and I will not give 
the formula. However, by formally solving ( [21b|) as 



B, 



PT 1 + (f 2 {Rn+l +Rn + Rn-l + Rn-2 + Rn-s) + {q' 2 ~ 1) £ R 



■J ' 



j<n— 4 



we see that the case of (|15|) is included for 7 — > 0. 

On the other hand, from the orthogonality condition, performing the calculation (|17D one 
derives relations between the coefficients A n , B n , C n . For A n and C n one finds again ( fflj ) 
for the special choice a = respectively 7 = (1 — q^ 1 )q A (2) q -i(4) q -iab. Furthermore, for B n 
we obtain 



Bn + (<? — 1)? 2n+A lRn-2Rn-?,Rn-AB n -2 

= g 2 (4) g - 1 6 7 - 1 + q- 2 (R n+1 + R n + # n _i + Rn-2) + q 2 ~ n R 



n—3 5 



(22) 
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which turns out to be consistent with (|2l), and in addition we have the following equation 



n Wfl 



(2) q -ia + (4) q -ib(R n+1 + R n ) - 7 g" n J R n _i J R n _ 3 



+ 1Q~ A [Rn+l (Rn+2 + Rn+1 + Rn) + Rn (Rn+1 + R n + Rn-l)] 



+ lBr, 



B 



n+2 



q 2 (4:) q -ib-f 1 - q~ 2 {R n+3 + R n+2 + R 



n+l 



Rr, 



+ (1 — 5)7 9 " R n Rn-lRn-2Rn-2,Rn- \ 



(23) 



which is the actual g-deformed string equation (at least for the lower order terms) , identifying 
(3 = q 2 (A) q -\b. Also eq. (|23|) is consistent with fl2~l|), which can be checked explicitely by 
tedious algebra. Thus, the orthogonality is consistent with the compatibility condition, 
which in turn implies that the g-orthogonal polynomials are indeed solutions of the system 
(PI). 



5 Discrete-Time Flows 



The dependence on the 'coupling constants' a,b, . . . as time-flow parameters in ( ]T5| ) is sig- 
nificant, and will lead to a discrete-time version of the Volterra system. In fact for the 
polynomials P n , we can derive from (|l6|) the g-difference time-evolutions of the form 



2D P 

q*DbPn 



QnP'n—2 1 

U n P n -2 + VnPn-A 



(24a) 
(24b) 



where the coefficients Q n , U n and V n are determined by the orthogonality condition. If 
we would have included additional higher-order exponents in the weight function of (|I"6|) 
this would lead to a discrete-time Volterra hierarchy. Let us illustrate this by working out 
explicitely the lowest order time-flow, namely in terms of the variable a. From ( j24a|) together 
with ( |lla| ) we obtain the relations 



q 2D a R n — Q n — Q n +l 
q 2 TaRn)Qn—l QnRn—2 1 



from which one can derive that 

Qn cRnRn—l ) Rn 



Rn 



{1-q )aQ n+l 



constant . 



(25a) 
(25b) 



(26) 



This then leads immediately to the following exact discrete-time g-deformation of the Volterra 
equation 



2D a R r , 



RnRn- 



q q 2 T a (R n+ iR n ) 



(27) 



It is slightly surprising that the exact integrable discretization of the Volterra systems 
turns out to take a form which resembles closely the original continuous-time equation, 
provided one considers the proper variables R n . Furthermore, solving ( 25b ) by taking 
Q n = ( q 2T a h n )/h n -2 , recalling that R n = h n /h n -i , we are led to the following equa- 
tion for h„ 



h„-i(q2T a h n ) = h n + (1 - q 2 )a( q 2T a h n +i) K-i + (1 - q 2 )a( q 2T a h n ) 



(28) 
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where taking c = 1 is consistent with the orthogonality. On the other hand, from the 
orthogonality condition flTBI ) one can derive by a calculation similar to flT7p 



2^T a h Ti 



q ,D a h n = q ,T a h n+1 + ] + (1 - g W > ( 29 ) 



which can be shown to be consistent with (|28|) . Let me mention that eq. (^) is a special 
subcase of a universal lattice equation derived some years ago in [^1 . 



Of course, the above derivation can be extended to the higher-order time-variables, e.g. 
starting from ( |24b| ). In general this will lead to more complicated systems and we will abstain 
in this note from presenting their derivation. What is important to note, however, is that for 
all the coefficients a, b, ... in the weight function (which in certain contexts are interpreted as 
coupling constants), the g-deformation in terms of the spectral parameter x leads necessarily 
also to a g-deformation in terms op these higher-order 'time'-variables which in turn leads 
to integrable discrete-time systems like the discrete Volterra system of (pTl). 



6 Discussion 

I have shown that it is fairly straightforward to obtain a g-analogue of the discrete Painleve 
I equation by straightforwardly g-deforming the continuous isomonodromy problem for dis- 
crete PI. In this way one obtains eq. ( p!4|) which is associated with probably the simplest 
isomonodromic deformation problem of a linear g-difference equation. Having this system 
at our disposal, we can now seriously embark on the more serious problem of investigating 
the full asymptotics for the corresponding transcendents. 

Furthermore, I have put forward a connection with g-orthogonal polynomials. For this 
one needs to consider a slightly more complicated system including 5th order terms in the 
"string equation". The connection with orthogonal polynomials is interesting, because on 
the one hand it establishes what is, in fact, a special similarity solution of the discrete- 
time Volterra hierarchy. On the other hand, it suggests a discretization of the Hermitean 
one-matrix model, which -on the continuum level- is well known to be related to orthogonal 
polynomials of Painleve type. It is suggestive, therefore, to try to trace back the g-orthogonal 
polynomials to the corresponding matrix models. If one were naive, one would be tempted 
to write down an expression of the form 

Z N (a,b) = [ [d q H\e-^ H \- u l H \ (30) 

JNxN h h 

for the g-deformation of the partition function of the Hermitean one-matrix model, in which 
[d q H] denotes a proper g-analogue of the Haar measure. However, eq. ( p0[ ) cannot be correct, 
because the g-exponents do not decompose in a natural way for linear combinations in their 
arguments, nor do the Jackson integrals allow for arbitrary changes of variables. What one 
needs is an interpretation of an expression of the form of (^) which should lead -after 
integration over the "angle variables" associated with an Hermitean matrix H to a multiple 
Jackson integral over the eigenvalues of H of the form 

N 

II (A* - A,) 2 • (31) 

i<j=l 



1\ 

n 



d q Xi 



-oA? - 

,-2 e„ 



-4 
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It is this expression for a 'partition function' that has a direct connection with the qPI 
equation of section 4. However, how to arrive at this expression starting from a (discrete) 
matrix integral needs some further study which is beyond the scope of this paper. 
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